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16 Abstract 

The Video Inertial Pointing (VIP) System is being developed to satisfy 
the acquisition and pointing requirements of astronomical telescopes. A 
unique feature of the system is the use of a single sensor to provide infor- 
mation for the generation of three axis pointing error signals and for a 
cathode ray tube (CRT) display of the star field. The pointing error signals 
are used to update the telescope's gyro stabilization and the CRT display is 
used by an operator to facilitate target acquisition and to aid in manual 
positioning of the telescope optical axis. 

The VIP system development includes analysis, simulation and the devel- 
opment of prototype system hardware. An early model of the system using a 
low light level vidicon was built and flown on a balloon-borne infrared 
telescope and is briefly described. The advanced VIP system, currently 
under development, can accept inputs from either a vidicon type sensor or 
from a state of the art charge coupled device (CCD) sensor. The advanced 
system hardware is described and an analysis of the multi-star tracking 
and three axis error signal generation, along with an analysis and design 
of the gyro update filter, are presented. Results of a hybrid simulation 
are described in which the advanced VIP system hardware is driven by a 
digital simulation of the star field/CCD sensor and an analog simulation 
of the telescope and gyro stabilization dynamics. 
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A COMPUTER-AIDED' TELESCOPE POINTING SYSTEM' UTILIZING A VIDEO STAR TRACKER 


James P. Murphy, Kenneth R. Lorell, and Charles D. Swift 
Ames Research Center ' 


INTRODUCTION 


Development of the Video Inertial Pointing (VIP) System was initiated at 
the Ames Research Center in 1973. The system concept and development has its 
origins in infrared (IR) astronomy investigations currently underway and the 
concept includes features which anticipate the acquisition and pointing 
requirements for astronomical telescopes of the Space Shuttle era. 

The initial version of the VIP system was developed to satisfy the 
acquisition and pointing requirements of a balloon-borne IR telescope. This 
system used a conventional SIT vidicon boresighted to the telescope to pro- 
vide information for a cathode ray tube (CRT) display at the ground station 
operators console and information to update the telescope’s gyro stabiliza- 
tion. An operator at the control console completed acquisition manually using 
a joystick and the CRT display. He then was able to select a single guide 
star for use in automatic on or off axis tracking. For this paper, on axis 
tracking refers to tracking a visible target on the common optical axis of 
the video sensor and telescope. Off axis tracking refers to tracking a 
visible target in the field of view (FOV) of the video sensor but not on the 
common optical axis. This system was used for several balloon flights and 
proved to be very useful. 

Because of the value of the early system, it became apparent that the 
concept had application to the more advanced airborne and Space Shuttle tele- 
scopes, such as the one meter class cryogenically cooled IR telescope being 
considered as a Shuttle payload (ref. 1>. Consequently, a program was ini- 
tiated to develop the VIP system to satisfy the requirements of these more 
advanced applications. Requirements of these telescopes include precise 
three axis pointing, manual and automatic modes of operation, and an interface 
with more’ complex and precise gyro stabilization systems. 

Key aspects of the advanced VIP system include use of a state-of-the-art 
charge coupled device (CCD) video sensor >which provides the capability for 
precision pointing; multi-star processing algorithms to' generate the three 
axis pointing error signals and a. gyro update filter which combines the video 
and gyro sensor data. The analysis and development of the multi-star proces- 
sing and gyro- update filter are presented along with some system tradeoffs 
that must be considered for design of an operational system. A prototype 
of the advanced system has been constructed and tests of this hardware with a 
hybrid simulation of the CCD sensor and telescope dynamics verify the system 
design and performance. Laboratory tests with a prototype CCD sensor have 
been conducted and the results are used in the analysis and hybrid simulation 
(ref. 2). Although not reported in this paper, the prototype CCD sensor will 



be used with the VIP system hardware on a ground based gimbal and gyro stabil- 
ization system for complete system hardware tests. 

It should be noted that the VIP system does not Include design or devel- 
opment of the gimbal and gyro stabilization. VIP is Intended to interface 
with this hardware and consequently, the design aspects of the gimbal and 
gyro stabilization are not discussed in this report. Similarly, for tele- 
scopes that require ultra precise pointing accuracy (below 1 arc sec) or that 
have significant movement of the optics with respect to the outer telescope 
structure, an internal fine guidance sensor is required. Although this 
report only treats a system using an external sensor boresighted to the tele- 
scope, the techniques can be applied to a system using a video sensor for 
fine guidance. 

The authors gratefully acknowledge the assistance of Messrs. W. F. 

Barrows of NASA/AIC, J. K. Lee of the Informatics Co., .and A. J. Throckmorton 
of Stanford University in developing the VIP hybrid simulation and the many 
hours spent making it operational. 


CONCEPT OF VIDEO INERTIAL POINTING 


Various approaches have been used in the past to provide the target 
acquisition and stabilized pointing for astronomical telescopes (ref. 3). 

The rotation of the earth and its stability are used by ground-based tele- 
scopes with a knowledge of the target location in the celestial sphere. A 
"finder" scope allows the astronomer to observe the star scene and to make 
final adjustments to the telescope angles and drive rates. For nonground- 
based telescopes where operation is remote, gyroscopes are used as the primary 
inertial reference or stabilization and image disector type star trackers are 
employed to correct for inherent gyro drift. Typically, several of these 
star trackers are used (at least two to provide three axis pointing informa- 
tion) but they cannpt provide simultaneous pointing error information and 
Information for a star field display. Consequently, conventional systems for 
balloon, air- and space-borne telescopes do not provide the same interaction 
and control available to the ground-based astronomer. 

The Video Inertial Pointing System takes advantage of recent improvements 
in video sensor and electronics technology to provide a telescope pointing 
system with fewer sensors than the conventional system and one that allows a 
more direct interaction with the control through a star field display at a 
control console. The system uses both gyro-stabilizatio.n and a video type 
star tracker. Primary stabilization is provided by conventional inertial 
techniques. The telescope remains stabilized, within performance tolerances 
of the gyros and servo electronics, in the inertial coordinate system. Fine 
pointing and control are made possible using star field- information from the 
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video sensor attached (usually boresighted) to the main telescope. The guide 
star locations from the video sensor are compared with the previously-stored 
locations of these same stars. Differences between the current and desired 
locations of these stars causes the system to generate error signals that are 
used to update the gyro outputs. The telescope (and video sensor) thus move 
to a vector so that the desired and current locations of the guide stars 
coincide. Since the system can be made to come to a stabilized null with the 
guide stars at any desired location in the field of view, it can be used to 
track any point in the field of view even if there is no visible star at that 
point (i.e. , off-axis pointing). In addition, the system has information 
from more than one star so that it can provide three axis pointing control 
using only a single video sensor. 

Another feature of the VIP system is the star field presentation avail- 
able to the observer/operator . Since the star sensor is of the video type, 
it is only necessary to present the output on a CRT display. In contrast, 
the conventional star tracker must resort to electronic or mechanical scan- 
ning in order to build up comparable star field information. It is not pos- 
sible to provide simultaneous pointing error signals and star field display 
information. In the systems described in this report, direct high bandwidth 
video is not used to drive the CRT display. Rather the video sensor data is 
processed to allow more elaborate display options (in the case of the advanced 
system) and to reduce the bandwidth for acceptable transmission (in the case 
of the first generation system discussed in the next section) . 


FIRST GENERATION SYSTEM 


A first generation model of the Video Inertial Pointing System has been 
completed and integrated into a balloon-borne telescope system for testing and 
evaluation. Figure 1 is a block diagram showing the key elements of the 
system. 

The telescope is gyrostabiiized in two control axes perpendicular to the 
telescope line of sight, .elevation and cross elevation. Mounted on the tele- 
scope and boresighted to it is a conventionally scanned silicon intensified 
target (SIT) vidicon camera. Its 75 mm f/1.4 lens results in a field of view 
of approximately 7° x 10° and sufficient sensitivity to track +6th magnitude 
stars. Star field video and synchronization pulses from this camera are fed 
to on-board processing electronics consisting of a video conditioner, gate 
generator and coincidence detector. The video conditioner circuit generates 
a timing pulse for each star signal above a commendable threshold, and in 
addition, generates x, y coordinate signals and Intensity signals for up to 
ten stars. These values are telemetered to the ground station CRT display. 

The gate generator uses the TV synch pulse to generate four timing gates which 
can be moved (by telemetered command) relative to the synch pulse so that they 
can be positioned to surround any location in the video format . The duration 
of the gates can be varied to include as few as 5 or up to 64 video scan lines. 
The square dead area surrounded can be varied from 2 to 20 video scan lines. 
Star timing pulses and 'gate signals are fed to the coincidence circuit. If a 
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timing pulse is coincident with one of the gate signals then an appropriate 
output signal is coupled to the gyro servo system which moves the telescope 
in the correct direction to move the star's timing pulse into the dead area 
in the center of the timing gate pattern. 

Control of the tracking system is by commands telemetered from the ground 
station. The operator or observer views the star field on the CRT display. 
Intensified cross hairs are visible on the display as well. The location and 
size of these cross hairs indicates the location and size of the timing gates 
in the video format. By manipulating a joystick, the cross hairs can be 
moved to cause a coincidence with the desired star. The on-board system will 
then move the telescope to put the star in the center of the pattern. Fur- 
ther movement of- the cross hairs is then possible to put the tracked star at 
whatever location is required in the field of view of the TV camera. 

Since this first system has only two axes of stabilization, periodic 
adjustments to the tracking are required to correct for diurnal motion. These 
adjustments are only necessary when the tracked star is not at the center of 
the field of view (off axis tracking) . During ground-based tests and two 
flight tests this system performed very well. Although the tracking accuracy 
depends on whether the guide star is on axis or off axis, the object of 
interest was within one arc minute of the optical axis of the telescope for 
extended periods of time during typical operation. 


ADVANCED SYSTEM DESCRIPTION 


The advanced VIP system consists of three primary subsystems; a video 
sensor whose optical axis is coaligned with the telescope optical axis, an 
electronics package which includes a microprocessor and interface electronics, 
and a control console with CRT display from which the system is operated. 
Figure 2 depicts, the inter-relationship between the telescope and VIP systems. 
In this section a description of the VIP components will be given, followed 
by a discussion of system operation. 

The video star tracker utilized by VIP is a key element in the system. 

Two important features of VIP are dependent on the video sensor; the ability 
to provide three-axis control with only a single star tracker, and a CRT 
display of the telescope FOV for the instrument operator. The advanced VIP 
system is designed to accept inputs from two types of sensors, either a con- 
ventional low-light level camera utilizing an Intensified SIT (ISIT) vidicon 
tube, or a state-of-the-art CCD camera which is all solid state and incor- 
porates its own microprocessor. While the ISIT sensor is a proven device and 
embodies relatively little technological risk, it does require an additional 
interface to the VIP computer in order to provide the required digital input 
format for multiple star tracking. The CCD tracker, in addition to the bene- 
fits accruing to an all solid-state device, has remarkable flexibility and is 
ideally suited to the VIP system. Both types of sensors employ interpolation 
techniques in order to obtain improved resolution. The y-processor in the 
CCD sensor electronics contains an algorithm to do digital interpolation. 
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The ISIT camera, interface electronics also contains a y-processor in order to 
perform the interpolation and data formatting functions. A more detailed 
description of the CCD tracker and its operation is presented in a later 
section. 

The VIP electronics package has two main sections, a National Semiconduc- 
tor IMP-16 li-processor and a group of circuits which interface the p-processor 
to the star tracker, telescope control system, and VIP control console. 

Figure 3 is a block diagram of the VIP electronics. 

The p-processor plays a central role in the VIP system. In addition to 
converting star tracker data into three-axis error signals via the multistar 
processing algorithm, it performs all of the logic sequences necessary to 
initialize and track, it drives the CRT display, and interrogates VIP control 
console switches to perform operator-controlled functions. The VIP 
p-processor also communicates directly with the star-tracker p-processor in 
order to effect data transfer as well as control the integration time and 
threshold level of the tracker. Figure 4 is a flow chart of the program used 
by the p-processor. 

In general, the entire electronics subsystem is designed to minimize 
hardware and place as many logic, data manipulation, and sequencing functions 
in software as possible. Thus, the primary function of the circuitry sur- 
rounding the IMP-16 in figure 3 is to control Input/output. Inputs from the 
control panel are detected by two jump-condition lines fed by a multiplexer 
which decodes the device address. Star— tracker data is clocked directly into 
the 16-bit data bus. Outputs are handled by 8-, 12- , and 16-bit latches 
which store the data when their appropriate addresses are decoded. CRT dis- 
play requirements dictated that a special buffer be used to store video -output 
data. One 12-bit D/A is used to drive analog circuitry by storing its output 
in sample and hold circuits selected in the correct sequence by the address 
decoder. Additional circuits are used for switching, impedance matching, and 
filtering. 

The VIP control panel and CRT display are shown in figure 5. The opera- 
tor controls the position of the telescope LOS with the joystick and may 
observe' a variety of star-field .presentations on the CRT display. Switches 
controlling VIP operation, auto/manual select, fast/slow slew rate, and star 
field display options, surround the joystick controller. Threshold level and 
integration period for the star tracker are selected by means of two thumb- 
wheel switches on the operator’s left. 

Other functions, such as automatic scan, are controlled from switches 
on the operator's right. Although the prototype equipment has an IMP-16 
which is programmable with paper tape, an operational version would use a 
p-processor with PROM program storage and would be co-located with the rest 
of the electronics package. 

Operation of a VIP-controlled telescope begins with coarse acquisition 
of the target star field. This is accomplished by slewing the telescope to a 
precomputed orientation and making an initial alignment with the joystick 
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controller in the manual mode. The star field, as viewed by the video sensor, 
is displayed on the CRT and the operator orients the telescope by comparing 
the CRT picture with a star chart. Alternatively, a computer-stored star map 
of the desired viewing area may be presented. When the target has been iden- 
tified, the automatic mode is selected by pushing the "initialize" button 
located on the end of the joystick. This commands the VIP p-processor to 
store the x~y locations of all of the stars being tracked by the video 
'sensor. Each time the current star locations are transmitted to VIP by the 
star tracker, the pointing error in all three axes is calculated by comparing 
these current star positions with those stored in memory upon the "initialize" 
command. The calculated attitude error and the gyroscopically measured atti- 
tude error are then combined in the y-processor gyro update filter to provide 
the torque commands fed back to each of the three gyros.- 

The operator may continue to reposition the telescope optical axis during 
automatic operation. A display option may be selected in which only the star 
field in its desired location is presented. Joystick activation moves this 
pattern on the CRT screen to whatever new position is desired. The dynamics 
of telescope motion are not seen (since only commanded star field position is 
displayed) thereby enabling very precise positioning. Simultaneous selection 
of current telescope position for the display will indicate to the operator 
how close the telescope is pointing to the desired location. The CRT will 
show two star fields in this mode, a bright one for actual position and a 
dimmer field for the desired position. When the two fields are coincident, 
the telescope is on target. 

A digital indicator on the control console informs the operator of the 
number of stars currently being used for tracking. Another display option 
permits viewing only these stars , As new stars enter or depart the field of 
view, VIP has the capability to display them without disturbing the tracking 
function, as well as continually keeping the operator informed of the number 
and location of tracked stars. If, because of slewing or tracker parameter 
changes, the number and/or position of tracked stars is unsatisfactory, the 
operator may reinitialize the system by pressing the "initialize" button 
again - 

VIP is designed to permit maximum operator interfacing with the auto- 
matic control system and provides him with several display options, indicators, 
and controls so that precise orientation of the telescope becomes a simple 
task. 


CCD SENSOR DESCRIPTION 


The advanced VIP system is designed to use either an ISIT vidicon camera 
or a sensor employing a charge coupled device (CCD) detector. Although con- 
ventional sensors do provide adequate performance, as demonstrated in the 
first generation system, there are a number of advantages, particularly for a 
digital system, that accrue from the use of a CCD sensor. Some of these will 
be discussed below, followed by a description of CCD sensor operation. 


6 



The primary advantage of a CCD sensor, compared with trackers utilizing 
vidicon or image dissector tubes, is the elimination of the scanning electron 
beam and the accuracy with which digital position information for multiple 
stars is generated. The absence of the drive circuitry, electron gun, control 
yoke, and other components associated with a scanning beam device is a con- 
siderable simplification in star tracker technology. Since the sequencing, 
control functions, and computation in VIP all rely on a digital computer, the 
inherent digital nature of a CCD sensor means the elimination of an analog- 
to-digital interface giving a considerable savings in system complexity. In 
addition, since all components are solid state, there are no high voltages, 
the operational life and reliability are enhanced, and the weight, size, and 
power all reduced. 

The face of the CCD detector is divided into a rectilinear pattern of 
discrete photosensitive elements, each connected to readout registers also 
located on the front surface. As photons strike the detector, the charge 
built up in each element is transferred to the readout registers after a 
selectable period of integration, the charge being directly proportional to 
the number of photons striking an element during the integration period. 

Since the location of each element is uniquely specified, it is possible to 
obtain' the precise coordinates, within the resolution of the photosensitive- 
element matrix, of any stars being imaged on the surface of the CCD. 

The current state of the art for CCD's is a matrix of about 400 elements 
square. Commercially available CCD's are about 200 elements on a side. Thus, 
even the most advanced units cannot obtain resolution better than 0.5 arc min 
for 3° FOV optics. However, it is possible through a process of Interpolation 
to Improve the resolution by an order of magnitude (ref. 2). 

The interpolation scheme requires that the output of a four by four 
matrix of elements be sampled and averaged each time the CCD is scanned. 

This array is centered on the element which first exceeds a pre-selected 
threshold during the initial scan of the device. Thus, for each star image, 
there is a corresponding four by four matrix for which the charge "center of 
mass" is computed. I,t is possible, therefore, to obtain computed star loca- 
tions which are not centered on an element. The effect of this procedure is 
to increase the resolution by artificially increasing the element density. 

The interpolation algorithm and the data formatting for use by the con- 
trol system are of a complex nature and require extensive data handling for 
their implementation. This is accomplished through the use of a microproces- 
sor built into the electronics of the star tracker. Raw data from the CCD is 
read at the end of each integration period by the microprocessor which then 
computes the locations of the ten brightest stars in the field of view. ’ Each 
star is assigned an identification number and the four pieces of data for each 
star (x, y, brightness, i.d.) are transferred to the VIP computer when the 
star tracker microprocessor completes the data for all of the stars in a frame. 
Logic within the microprocessor automatically takes care of stars entering or 
leaving the field of view so that while it is always possible to display all 
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stars in the field of view, the VIP controller only accepts data for valid 
stars. ^ Initializations of the star tracker may be accomplished at any time 
at the option of the operator via a pushbutton located on the VIP control 
stick. 

The integration period and threshold are also controlled by the star 
tracker microprocessor. Seven levels of both functions are selectable from 
thumbwheel switches on the VIP control panel. These levels are arranged so 
that each Increment of integration time changes the sensitivity by 0.25 
and each threshold increment changes the sensitivity by 1.0 Thus, if 

= KT^ (2.5118~“^) (1) 

describes the charge built up in a given element during an Integration period 
j, and we desire to have a +0.25 change for each increase T j , then 

Q(Tj,M^) = Q(T^_j_^,M^ +0.25) 

Tj+l (^.5^1®"^^^^°’^^^) = T^ (2.5118”^) 
or 

- 1.2588 Tj (2) 

For a given Tj , Qj^(Tj), the threshold below which a star will be rejected as 
being too dim, may be expressed as 


= (2.5118)“1q^(T^) (3) 

From equations (2) and (3) we have a method of computing the appropriate inte- 
gration times and threshold levels given the minimum exposure time, system 
gain constant K, and the maximum desired M^. For practical considerations, 
^^min ~ 0-250 sec and =8. It is important to note that the relation- 

ship between the actual integration times or threshold levels and the posi- 
tions of the thumbwheel switches on the control panel are arbitrary, the star 
tracker microprocessor need only use the switch position code to locate the 
desired value of integration time or threshold in a look-up table stored in 
memory. 

Table 1 lists the values of Qi(Tj) and the range of (brightest and 
dimmest stars) for any combination of threshold and integration time. K has 
been arbitrarily set to 10^. 


valid star is one that has met a number of tracking criteria, both in 
the star tracker and within the VIP software. Attitude information is gener- 
ated using only valid stars. 
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Figure 6 is a view of the star tracker cut-away to show the major com- 
ponents. A modular approach has been taken to the optical and mechanical 
design so that the tracker may be easily modified to fit a variety of mission 
requirements. A Peltier-effect cooler is included for CCD operation on the 
ground; operation of the cooler during high altitude or orbital flight is not 
required. 


VIP GYRO FILTER ANALYSIS 


A number of investigators have approached the problem of combining star 
tracker information and gyro output to improve inertial orientation. Farren- 
kopf (ref. 4) uses a precise gyro error model with pre— and post— update filter 
error analysis to construct a family of curves which predict the performance 
of systems of this type. Rauch (ref. 5) has analyzed the errors in a single- 
axis star-tracker/gyro system. He obtains the rms steady state pointing error 
as a function of the number and accuracy of star measurements and the accuracy 
of the gyro. Potter and Vander Velde (ref. 6) derive an analytical solution 
to the problem of combining gyro and star tracker data in an optimal fashion. 

In all of the work described above, the authors have investigated the 
interrelationship between high .bandwidth, high (relatively) noise gyros and 
low bandwidth, low noise star trackers when their outputs are combined for 
inertial stabilization. However, the problem of design, and in particular 
the interfacing of digital and continuous systems, has not been addressed. 

In this section the design of a simple filter for the VIP system is discussed - 
and the results of an analytical technique which utilizes computer-generated 
Z transforms is presented. 

Figure 7 is a block diagram illustrating the basic elements of a gyro- 
stabilized telescope controller to which a VIP system has been added. The 
telescope control loop is closed by coupling the output of a single axis rate 
integrating gyro to a torque motor through a conventional compensation net- 
work. The system has excellent transient response and, with integral control, 
has good performance even when a body-fixed disturbance torque is applied. 

The gyro output is updated by reference to the attitude error measured by a 
star tracker. 

The Separation Theorem Implies that the operation (i.e. , pole placement) 
of the update filter L(S) will not affect the plant controller as long as the 
filter model reflects the plant dynamics. Thus, if a continuous measurement 
of attitude error were provided, the design of the gyro update filter could 
proceed along the lines of classical control theory. For example, since the 
goal of the controller is to use 0q to force G = 0 and that of the filter 
is to use y to force 0q = 0, then 

Hc[dg(S)GD + 

JS2(1 + LGD) + HCGD(S + L) 


0(S) = - 


(4) 



( 5 ) 


- 

for a simple gyro approximation 


e(s) 

'eg(s) 


[dg(S)GD + 6gJjs2 
JS2(1 + LGD) + HCGD(S + L) 


GD . i d^<S) - -f so 



(JS2 + HC) (S + L) 


■J(“g^ + °Go^') 

(JS2 + HC) (S + L) 


thus equations (6) and (7) imply a simple filter of the type 


( 6 ) 


(7) 

^2 

L(S) = Ki + — 


will eliminate d and An alternate method of arriving at the same form 

is to mechanize an observer for 6g using 6g as a measurement. Then the 
system dynamics are 


®G = *^G. + ® 

d_ = 0 
G 

and for an observer 

«G - ^G + - ®G> * 

^G - - 8g> 


<S) 


(9) 


Equations (9) lead to the same form of filter as before. Although typical 
systems of the type shown in figure 7 have 10 or more states, the procedure of 
selecting Ki and K 2 is a straightforward task for a linear system. 


The use of a digital star tracker and digital computation, to determine 
the attitude error introduces a number of elements into the system which have 
a direct, and adverse, effect on the performance. These include sampling, 
time delays in both the star tracker and computer, mixing .of delayed and non- 
delayed data, and quantization. Thus, a major difficulty is introduced into 
the analysis and design of a system with these characteristics. The update 
filter couples into the stabilization loop dynamics and unstable operation of 
the two loops will result if the gains and K 2 are selected without pro- 
perly taking into account all of the dynamical effects present. 

Classical S-plane analysis is still possible for a system with sampling 
and time delays if some form of polynomial representation of such as 

the Fade approximant (ref, 7), is used. This technique is mathematically 
cumbersome, however, and provides relatively imprecise knowledge of the gain- 
response interaction. An alternative is the Z-transform which is ideally 
suited for the analysis of problems of this type, but has the disadvantage of 
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being difficult to calculate, especially for systems with a large number of 
states. This problem has been circumvented through the use of computer pro- 
grams (ref. 8) which calculate the Z-transform, and find the roots of the 
resulting characteristic equation. 


A two-step process is used to obtain the $ and F matrices of the com- 
plete system. First the continuous portion of the system, without time delays 
or digital filter, is transformed. Augmented state transition and control 
distribution matrices may now be formed by adding the appropriate Z-transforms 
for the time delays and digital portion of the control (filter). Thus, for 
step one: 


X = Fx + Gu -3- X , , = + Fu 

— — — — n+1 — n — r 


and for step two: 
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X 
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XX 

xy 
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y 1 1 
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y 
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yx 

yyj 


■^n 


( 10 ) 


Z— plane plots of the eigenvalues of equations (10) , analogous to closed loop 
root loci in the S-plane for a continuous system, permit the designer to 
select filter gains to optimize whatever criteria he desires. In addition, 
the changes in performance as a function of the time delays may be observed 
and corrected for. In> particular, T]^, the CCD integration time, has a wide 
dynamic range. 


Figure 8 is a state variable representation of the system shown in 
figure 7 with numerical values for a gyro stabilized telescope used at the 
Ames Research Center. This system, with appropriate values for the time 
delays and sampling intervals, has been analyzed. A filter similar in form 
to the one discussed for the case "of a linear system was implemented in the 
form of the difference equation 


J'Gn = >'G(n-l) ' Vcn-i)' ' »» 

L(Z), the Z-transform of the filter is 

(Z - V 

- >1 (Z - 1) »2) 

where 

\ (Ki + K2) , Kp ^ 

Figures 9 and 10 represent a series of typical Z-plane loci calculated 
for this system. The effects of gain changes and variations in the delay 
times are quite evident. The general effects of the time delays is to 
strongly increase the sensitivity to changes in and require much smaller 
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values of K 2 . A natural result is that speed of response of the entire 
system Is considerably lessened and transient behavior is degraded. 

A conclusion to be drawn from this analysis is that a more sophisticated 
filter mechanization which models the time delays is required in order to 
improve performance. While our prototype system utilizes the simple filter 
described above, investigations are currently under way into the analysis, 
design, and implementation of filters which incorporate a more precise plant 
model . 


VIP MULTI-STAR PROCESSING 


To determine the attitude error signals needed to point a telescope, at 
least two guide stars must be used. If a single guide star were used pointing 
information for two axes perpendicular to the telescope’s line of sight could 
be .generated. Gyro information, subject to drift, would then be required to 
orient the telescope about its optical axis. Since the video sensor used with 
the VIP system provides position information for up to 10 guide stars, the 
complete three axis pointing errors needed to point the telescope, as well as 
error signal averaging used to eliminate the effects of video sensor noise 
and inaccuracies, can be generated. 

Appendix A contains a derivation of the basic equations relating the 
desired position of the guide stars to the actual measured position and the 
resulting three axis pointing errors. A small angle approximation was used 
for this derivation since stars seen by the video sensor are within a small 
angular distance of the telescope optical axis and a more sophisticated deri- 
vation using direction cosines or quaternions would result in equations that 
require an excessive amount of computer time to solve. 

With the basic equations derived in appendix A, several methods were 
considered for deriving the final or best estimates of the three axis pointing 
errors. One approach considered was to use all the star data with a least- 
mean-squares approach assuming some error statistics associated with each star 
position measurement based on position in the CCD array and star magnitude. 

This method was compared to simply selecting star pairs, and solving the 
appropriate algebraic equations for the pointing errors. The least mean 
square approach resulted in more complex equations and, further, offered 
little benefit unless there was a great difference in the error statistics 
associated with the star measurements. 

Consequently the approach taken for the multi-star processing is to 
select star pairs and the appropriate algebraic equations and then determine 
the three axis pointing errors. The resulting pointing errors for each star 
pair are then averaged together to obtain a best estimate. 

Table 2 contains a summary of the algebraic equations for the pointing 
errors derived in appendix A. For each star pair there are two independent 
equations for roll and three independent equations for pitch and yaw. Assuming 
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identically distributed, independent statistics for the measurements of star 
position, the variances of the attitude error equations are also included in 
table 2. Each of the algebraic equations has singularities as well and a 
selection of the appropriate equation must therefore be based on the relative 
location of the two, stars. Table 3 summarizes the selection process. The 
only restriction occurs when two stars are in the same quadrant and the singu- 
larities in the equations and resultant amplification of error indicates that 
this case should be rejected as a valid star pair. Of course there is an 
added restriction for the selection of star pairs to calculate the roll point- 
ing error. Due to the geometry of using stars relatively close to the tele- 
scope line of sight, the roll pointing error is much larger than that for 
pitch and yaw and star pairs for roll should be selected as far from the 
origin within the field of view as possible. 

The ultimate performance of the multi-star processing depends, of course, 
on the performance of the video sensor and the number of star pairs used. 
Performance of the CCD sensor will include two effects: accuracy and quanti- 

zation of the interpolation calculations, and noise associated with random 
fluctuations of the levels in the individual CCD wells . Design of the CCD 
sensor ensures that the quantization effect in the interpolation calculations 
will be small and current test data shows that the interpolation calculations 
accuracy is approximately ±10 percent of the CCD well size. Noise associated 
with the CCD sensor star position measurements is also small compared to the 
interpolation accuracy and will not be considered in this section; however, 
it is important for design of the gyro filter. 

Since the interpolation error is the primary error source from the CCD 
sensor, a measure of the accuracy (neglecting misalignments) of the VIP system 
can be closely approximated as follows. If the interpolation error is assumed 
to be uniformly distributed then 

or Ay^^l = ~ = (13) 


where a is ten percent of the CCD well dimension in arc sec. If three guide 
star pairs, allowing use of equations (A5a) from table 2, are assumed. 


and 
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then 
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E(§^2) 



Similarly, for roll, if three guide star pairs with the stars at least I” 
from the center, are assumed. 


and 


(y - “ y -)^ 

'■'oi ^ 03 ' 


and 


1 , MO 

(x . - X .)2 ^ IT 
01 oj 


2) 180 X a ,,. ^ 

^ R ' - IT 


(15) 


and 


E(0j^2) 


< 180 X 2 X 
3it 


Table 4 summarizes the 1 a null pitch and- yaw pointing accuracies for both 
the 190x244 and 400x400 CCD array sizes. The poorer accuracy for roll is 
clearly shown and is a result of the geometry as explained earlier. While 
the roll numbers appear large, it must be remembered that this is accuracy 
and is not a measure of the roll stability due to quantization effects in the 
CCD sensor. The roll quantization would be much smaller. Table 4 clearly 
shows the effects and benefits of the multi— star processing. 


VIP HYBRID SIMULATION 


Operation of the complete VIP System electronics has been thoroughly 
tested utilizing a hybrid computer to simulate the functions of the CCD star 
tracker, the telescope dynamics, and the gyro-controlled stabilization system. 
Figure 11 is a block diagram representing the interrelationship of the major 
components of the simulation. 

An HP-2114 digital computer is employed to generate the star field, per- 
form the logic functions required to select stars that would be detected by 
the CCD star tracker, and format and output the data in a manner identical to 
the Intel 8080 microprocessor in the star tracker electronics. The telescope 
control system and dynamics are programmed on an EAI TR-48 analog computer. 

The body pitch, yaw, and roll rates are sampled directly from the analog com- 
puter by the digital computer which uses them to determine the apparent -motion 
of the star field. Simulated gyro outputs from the analog computer are com- 
pared to the computed pitch, yaw, and roll error angles calculated by VIP from 
the star tracker data. The results of this comparison is used to drive an 
update filter which applies the appropriate signals to the gyro torquers. 

Thus, the 2114 and TR-48 are interconnected and each has two-way communication 
with the VIP electronics package as well. 


14 



The following section provides a brief description of the CCD star 
tracker simulation on the HP-2114 computer. Figure 12 illustrates an 
inert ially- fixed coordinate system whose origin is at the telescope center of 
mass. Stars are located on an imaginary celestial sphere of unity radius 
centered about this coordinate system. The right ascension and declination 
for any star may be converted to the corresponding position vector by use of 
equations (16) below: 


X = cos 0 sin ^ 1 

o I 

yg = sin 0 > (16). 

z = cos 0 cos (j> / 

O 

The right ascension, declination, and magnitude of up to 30 stars are 
stored in the HP-2114 memory. The field has been limited to 16° square and 
several test patterns as well as constellations taken from star charts have 
been programmed. Figure 13 shows a typical pattern. 

The location of a star in the body-fixed coordinate system of the 
teles cope- star tracker can be calculated using the direction cosine matrix 
which relates the body and inertial coordinate systems. The elements of the 
direction cosine matrix are updated by the HP 2114 after each sample of the 
analog states representing the actual (as compared with the measured) body 
rates . Thus 


Xb = 


(17) 


Since the star tracker is only concerned with the x and y coordinates, 
the computer actually calculates 


(18) 


^b' 
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tests first each Xg, followed by the yg's of the stars that remain. So 
that if both 

-^F - ^ ^F 

1 

-yp ^ yg. ^ yp 

X 


( 19 ) 


the coordinates of that star are saved. 
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The integration time and threshold are selectable functions from the 
VIP control panel. During each data cycle, the status of these functions is 
sampled and the following algorithm is used to determine the of the dim- 

mest star capable of being -detected by the tracker; 

^ (6.75 - j) + 0.25i (20) 

where j is the threshold level and i is the integration- time switch posi- 
tion. Thus, stars are eliminated from the tracker field of view if they are 
dimmer than the maximum calculated in equation (20) , 

Two additional functions are required for the star tracker simulation; 
data formatting and time delay. These are both performed at the end of the 
computation cycle. Formatting is necessary because a) the coordinate system 
used by the star tracker has its origin in the lower right corner instead of 
the center, b) the data must be delivered to the VIP interface in the correct 
sequence and timing, and c) the data must be scaled from the Xg, yg, 
values used for calculation to the values the star tracker would actually 
output. The scaling simulates both the interpolation and quantization 
inherent in the CCD. A time delay in the real tracker results from the period 
of integration and the computation, necessary to do interpolation. The HP-2114 
logic automatically adjusts the time delay to the appropriate value according 
to the integration time requested by the VIP control panel switch. Figure 14 
is a flow chart for the HP-2114 program indicating the major steps required 
for the simulation. Figure 15 is a photograph of the experimental setr-up 
showing the major components. 

The results of tests using the hybrid simulation are shown in fig- 
ures 16(a) and (b) . The relatively slow response of the filter and the some- 
what under-damped nature of the commanded attitude response were predicted by 
the 2-plane analysis discussed in the previous section. The performance is 
more than adequate, however, and with the addition of the commanded star field 
display option, described in an earlier section, the operator may position the 
.optical axis without having to contend with the system dynamics. 


CONCLUSIONS 


Design and performance details of a first generation system and an 
advanced Video Inertial Pointing system capable of providing three— axis 
pointing error information with a single video sensor have been presented. 

The results indicate that the concept is a viable one. The use of a multi- 
star processing to generate three-axis error signals has also been shown to 
substantially improve star-position measurement statistics when a simple 
averaging algorithm is Included. An analysis of the problem of using delayed, 
quantized star tracker information to update the output of the gyro stabiliza- 
tion loop revealed the importance of careful modeling of the discrete 
portion of the system. A hybrid simulation in combination with work done at 
the Jet Propulsion Laboratory on the VIP CCD video sensor (see 2) has- 
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demonstrated the sucessful integration of state-of-the-art devices such as 
microprocessors and a solid-state video star tracker in an operating system. 

It is anticipated that the VIP system concept will he part of the Shuttle 
Infrared Telescope Facility as well as other systems used to point instruments 
on the Space Shuttle. 
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APPENDIX A 


DERIVATION OF ATTITUDE ERROR EQUATIONS 


Consider the small-angle approximation to the direction cosine matrix 
which transforms inertial to body-fixed coordinates through orthogonal pitch, 
yaw, and roll angles : 


1 

-e. 


9 -0 

R y 


(Al) 


0-0 1 

y p 

The position of a star in body coordinates after a small angle rotation is 




1 0„ -0 
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R y 


= 

-0_ 1 0 
R p 

z^ 


9-0 1 

B 


y P 


(A2) 


Since we're interested in the Xg, yg components only, we may expand them as 


-y“i 


■ -®R*I 


(A3) 


If x^, y^ is considered to be the initial location of the star, then the 
change in"^position, as seen by the star-tracker is 

Ax 4 x„ - X 
B 


i Yg - Y] 


(A4) 


Equations (A3) may be normalized by dividing through by Zj-, since 
» Xj,yj and the ^-j-'s are essentially equal for all stars. Defining K 
as the scaling factor for the' normalization and the conversion of the image 
through the electro-optical system of the star tracker, equations (A3) may be 
re-written as: 


KAx = 0„Ky - 0v 
R o i 

KAy = -e,,Kx + 

K O 


(A5) 
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Thus, for two stars, using the Ax equation 


and for the Ay equation 


Ax. - Ax. 

H = 1 J 

R (y . - y .) 

■'ox ^oj 
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03 ox 
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03 ox'^ 


An additional set of equations for 0p and By may be obtained by using 
the 0R derived from the Ay portion of (A5) in the Ax equation and vice 
versa. They are v 


K^fy -(^y- “ Ay.) + Ax.(x . - X .)1 
_ Pox -X ^3’' x^ ox 03 J 


(x . - X .) 
•03 ox 
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L ox X 3^ -X •'ox 03 J 
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Using equations (A6) and (A7) for a star field of n stars, we may 
improve our calculation of 0p^Y,R the results. Thus 
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We may also wish to utilize some form of weighted averaging based on, for 
example, the relative brightness of the star, the separation of a star pair, 
or the distance of the star from the telescope line of sight (for calcu- 

lations) . Given some weighting factor W.. = f(x.,y.,x . ,y .,x.,y.,x 
y ,,M ,M ) then equations (A 9 ) become i x . ox j 3 03 


03 V^ Vj 
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n n 
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TABLE 1.- STAR TRACKER DYNAMIC RANGE FOR SEVEN THRESHOLD LEVELS AND SEVEN 


INTEGRATION TIMES 


Thresh 

\s^vel 

sec 

TH 1 

TH 2 

TH 3 

TH 4 

TH 5 
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TH 7 
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Q > 0.40 
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TABLE 2.- EXPECTED VALUE OF ATTITUDE ERROR FOR EACH OF THE FIVE TYPES OF ATTITUDE ERROR EQUATIONS 



Attitude error equation 

E{Attitude error^}* 

Singularities 
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TABLE 3.- STAR-PAIR SELECTION CRITERIA 


Star-pair 

configuration 


Roll 

equation 

Pitch 

equation 

Yaw 

. equation 

Not selected due 
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TABLE 4.- 1-cT ROLL, PITCH, YAW POINTING ACCURACY 
(IN ARC SEC) WHILE TRACKING 4 STARS 

















(a) Airborne portion. ' 

Figure 1,- First generation VIP -system. 












GROUND STATION 



(b) Ground station. 
Figure 1.- Concluded. 
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Figure 2.- Advanced VIP system. 


















SYSTEM UTILITY DRIVERS 


(a) HP-2114 Star f ield/tracker simulation program structure. 
Figure 4.- VIP microprocessor software organization and flow chart. 











(b) VIP software flow chart. 


Figure 4.- Continued 














(c) VIP software flow chart 
Figure 4.- Continued. 
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(d) VIP software flow chart. 


Figure 4.- Continued 













(e) VIP software flow chart 
Figure 4.- Concluded. 
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Figure 6.- CCD star tracker. 
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Figure 7.- Block diagram of gyrostabilized telescope with VIP system. 












Figure 8.- State variable representation of gyrostabilized telescope with 

numerical values for AIROScope system. 
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Figure 11,- Hybrid simulation block diagram. 







Figure 


12.- HP-2114 star 


field coordinate system. 
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Figure 13.- Constellation Orion used in hybrid simulation 





(a) HP-2114 software organization 
Figure 14.- HP-2114 software organization and flow chart 


























(b) HP-2114 software flow chart 
Figure 14.- Continued, 
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(c) HP-2114 software flow chart 
Figure 14.- Continued. 
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(d) HP-2114 software flow chart. 
Figure 14.- Continued. 
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(e) HP-2114 software flow chart. 


Figure 14.- Concluded 
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Figure 15.- Laboratory set-up for VIP hybrid simulation. 
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Figure 16.- VIP system performance from hybrid simulation. Sample time = 0.25 sec, delay time = 0.25 sec, 
field of view = 2° x 2°, tracked stars =3. A. 10° /hr step drift input to roll and pitch gyros. 
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(b) Pitch axis filter response, 
'Figure 16.- Concluded. 


